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Abstract 



Natural ecosystems are characterized by striking diversity of form 
and functions and yet exhibit deep symmetries emerging across scales 
of space, time and organizational complexity. Species-area relation- 
ships and species-abundance distributions are examples of emerging 
patterns irrespective of the details of the underlying ecosystem func- 
tions. Here we present empirical and theoretical evidence for a new 
macroecological pattern related to the distributions of local species 
persistence times, defined as the timespans between local coloniza- 
tions and extinctions in a given geographic region. Empirical distribu- 
tions pertaining to two different taxa, breeding birds and herbaceous 
plants, analyzed in a new framework that accounts for the finiteness 
of the observational period, exhibit power-law scaling limited by a 
cut-off determined by the rate of emergence of new species. In spite 
of the differences between taxa and spatial scales of analysis, the scal- 
ing exponents are statistically indistinguishable from each other and 
significantly different from those predicted by existing models. We 
theoretically investigate how the scaling features depend on the struc- 
ture of the spatial interaction network and show that the empirical 
scaling exponents are reproduced once a two-dimensional isotropic 
texture is used, regardless of the details of the ecological interactions. 
The framework developed here also allows to link the cut-off timescale 
with the spatial scale of analysis, and the persistence-time distribution 
to the species-area relationship. We conclude that the inherent coher- 
ence obtained between spatial and temporal macroecological patterns 
points at a seemingly general feature of the dynamical evolution of 
ecosystems. 
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1 Introduction 



Understanding local extinction processes has gained urgency as the number 
of threatened species increases throughout the world because of factors like 
habitat destruction or climate change [U El El lU E] > but a synthesis of the- 
ory and empirical evidence accounting for the relevant ecological dynamics 
is lacking. In this context, we address here the study of persistence times of 
trophically equivalent co-occurring species in relation to the spatial scale of 
observation. The persistence time r of a species within a geographic region 
is defined as the time incurred between its emergence and its local extinc- 
tion (see P, [7] and Fig. [I]). At a local scale, persistence times are largely 
controlled by ecological processes operating at short timescales (e.g. pop- 
ulation dynamics, dispersal, immigration, contraction/expansion of species 
geographic ranges) as local extinctions are dynamically balanced by coloniza- 
tions [HI IH]- At a global scale, originations and extinctions are controlled by 
mechanisms acting on macroevolutionary timescales. 

From a theoretical viewpoint, the simplest baseline model for population 
dynamics is a random walk without drift, according to which the abundance 
of a species in a geographic region has the same probability of increasing or 
decreasing by one individual at every time step. According to this scheme, 
local extinction is equivalent to a random walker's first passage to zero, and 
thus the resulting persistence-time distribution has a power-law decay with 
exponent 3/2 [TO] . 

A more realistic description can be achieved by accounting for basic 
ecological processes like birth, death, migration and speciation in neutral 
[TTl [T2"l [T31 H3] mean field schemes, as follows. Consider a community of N 
individuals belonging to different species. At every time step a randomly 
selected individual dies and space or resources are freed up for colonization. 
With probability v the site is taken by an individual of a species not currently 
present in the system; v is equivalent to a per-birth diversification rate and 
it accounts for both speciation and immigration from surrounding communi- 
ties. With the residual probability 1 — u the died individual is replaced by one 
offspring of an individual randomly sampled within the community [151 IB] • 
As such the probability of colonization by a species depends solely on its rel- 
ative abundance in the community. The asymptotic behavior of the resulting 
persistence-time distribution (i.e. p T {t)) exhibits a power-law scaling limited 
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by an exponential cut-off: 

p T (t) oc t-"e-", (1) 

with exponent a = 2 [7J. In Eq. [IJ time is expressed in generation time 
units [TT], i.e. it has been rescaled such that the birth rate is equal to one. 
Notably, in the mean field scheme the probability distribution p T (t) depends 
solely on the diversification rate which accounts for speciation and migration 
processes and imposes a characteristic timescale 1/v for local extinctions. 
While per-birth speciation rates are not expected to vary with the spatial 
scale of analysis, per-birth immigration rates are argued to decrease as the 
spatial scale increases. In fact, the possible sources of migration (chiefly 
dependent on the geometrical properties of the boundary and the nature 
of dispersal processes) are argued to scale sub-linearly with the community 
size [I7J , which in turn is typically linearly proportional to geographic area [5J 
12]. As continental scales are approached, migration processes (almost) vanish 
and the diversification rate ultimately reflects only the speciation rate. 

From an empirical viewpoint, species and genera persistence times de- 
ducted from fossil record data have been suggested to follow either power-law 
(with non-trivial exponents in the range 1.5 — 2 [151 [19] [20J) or exponential 
distributions [2HHH]. It has been argued, however, that data quality, in par- 
ticular for species, precludes a critical assessment [TJ. Also, local analyses of 
species persistence over ecological timescales suggest power-law distributions 
with non trivial exponents [6]. 

In what follows we shall provide evidence for power-law behavior, either 
empirically or from a broad spectrum of theoretical derivations. Implica- 
tions on emerging macroecological patterns will be examined, with special 
attention to possible biogeographical effects. 

2 Empirical Persistence-Times Distributions 

We empirically characterize species persistence-time distributions by analyz- 
ing two long-term datasets covering very different spatial scales: i) a 41-year 
survey of North American breeding birds [22J; and ii) a 38-year inventory of 
herbaceous plants from Kansas prairies [23] . 

The North American Breeding Bird Survey consists of a record of annual 
abundance of more than 700 species over the 1966-present period along more 
than 5000 observational routes. The spatial location of the routes analyzed 
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is shown in Fig. 1. We consider only routes with a latitude less than 50° 
because density of routes with a long surveyed period drastically decreases 
above the 50th parallel. Noting that in many regions the survey started 
only in 1968, we discard the first two years of observations in order to have 
simultaneous records for all the regions in the system. The spatial extent 
of the observational routes allows us to analyze species persistence times 
at different spatial scales. We consider 20 different scales of analysis with 
linearly increasing values of the square root of the sampled area starting from 
A = 10000 km 2 to A = 3.8 • 10 6 km 2 . We also analyze the whole system, 
which corresponds to an area of A = 7.8 • 10 6 km 2 . For every scale of analysis 
A we consider several overlapping square cells of area A inside the system. A 
three-dimensional presence-absence matrix P is thus built. Each element p s t c 
of the matrix is equal to 1 if species s is observed during year t in at least one 
of the observational routes comprised in cell c, otherwise p s t c = 0. For every 
scale of analysis we discard the cells that (i) do not have a continuous record 
for the whole period (41 years) or (ii) have more than 5% of their area falling 
outside the system. For every cell and every species we measure persistence 
times from presence-absence time series derived from the second dimension of 
matrix P. Persistence time is defined as the length of a contiguous sequence 
of ones in the time series. For every scale of analysis we consider all the 
measured persistence times regardless of the species they belong to and the 
cell where they are measured. The effect of possible imperfect detection of 
species [24J on measured persistence times has also been investigated (see 
Supporting Information (SI)). 

The herbaceous plant dataset (23] comprises a series of 51 quadrats of 1 
m 2 from mixed Kansas grass prairies where all individual plants were mapped 
every year from 1932 to 1972. In order to meet the data quality standard 
required for our analysis as discussed above for the breeding bird data, we 
discard 10 quadrats and the first three years of observations. Due to the 
limited number of observational plots in the herbaceous plant dataset we 
limit our analysis to quadrat spatial scale A = 1 m 2 . Analogously to the 
previous case, we reconstruct the matrix P from presence-absence data for 
every species, year and quadrat. 

Note that, when dealing with empirical survey data, the effect of the 
finiteness of the observational time window on the measured species per- 
sistence times must be properly taken into account. To this end, we have 
developed new tools to extend the inference of persistence-time distributions 
for periods longer than the observational window. In particular, we ana- 
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lytically derive, given the persistence-time probability density function, the 
distribution of two additional variables that can actually be measured from 
empirical data: i) the persistence times r' of species that emerge and go lo- 
cally extinct within the observed time window AT W ; and and ii) the variable 
t" that comprises r 1 and all the portions of species persistence times that 
are partially seen inside the observational time window but start or/and end 
outside (Fig. |2]\ and Materials and Methods). The fmiteness of the time 
window imposes a cut-off to p T '(t). On the contrary p T »(t) has an atom 
of probability in t = AT W corresponding to the fraction of species that are 
always present during the observational time. By matching analytical and 
observational distributions for p T >(t) and p T "(t), it is possible to infer the 
persistence-time distribution p T (t). The scaling exponent and the diversifi- 
cation rate for the herbaceous plant persistence-time distributions have been 
determined with a simultaneous nonlinear fit of observational and analytical 
p T '(t) and p T "(t). Confidence intervals are equal to the standard error of the 
fit. For breeding birds, we repeat the nonlinear fit for different spatial scales 
of analysis. The reported scaling exponent and the confidence interval have 
been obtained by averaging results across spatial scales. 

Remarkably, the persistence times of breeding birds at different spatial 
scales of analysis and of herbaceous plants prove to be best fitted by a power- 
law distribution with an exponent a = 1.83 ± 0.02 and a = 1.78 ± 0.08, re- 
spectively (Fig. [2j3,C) . It is important to note that both scaling coefficients 
derived empirically are significantly different from the predictions of the ex- 
isting baseline models discussed above (the random-walk persistence time 
yields an exact exponent a = 3/2; the mean-field model yields a = 2). 

3 Theoretical Persistence-Time Distributions 

In this section we provide a theoretical rationale for the universality of the 
scaling behavior of persistence-time distributions with respect to the topol- 
ogy of the interactions allowed by the environmental matrix. In particu- 
lar we provide evidence on how non-trivial exponents of the type observed 
empirically can be reproduced by simple theoretical models once dispersal 
limitation and the actual network of spatial connections are taken into ac- 
count. We have implemented the neutral game described above in regular 
one-, two- and three-dimensional lattices in which every site represents an 
individual [13 [IB] . We have also explored the patterns emerging from the 
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application of the model to dendritic structures mimicking riverine ecosys- 
tems where dispersal processes and ecological organization are constrained 
by the network structure. Indeed, many features of riverine ecosystems have 
been shown to be affected by the connectivity of river networks [251 126] . In 
particular, river geometry has been studied in relation to extinction risk |27j, 
migration processes [28], persistence of amphibian populations macroin- 
vertebrate dispersal [30] and freshwater fish biodiversity [TU [31] . For general 
calculations of the topological structure and metric properties relevant to 
dendritic ecological corridors, we employ the features of Optimal Channel 
Networks (OCNs) [32J. They hold fractal characteristics known to closely 
conform to the scaling of real networks [33]. Among the advantages of the 
use of OCNs, one recalls the possibility to fit one such construct into any 
assigned domain (e.g. a square, Fig. [3j, thus allowing exactly the same size 
and number of nodes of a two-dimensional lattice to be endowed with al- 
tered directionality of connections. To account for limited dispersal effects, 
we allow only the offsprings of the nearest neighbors of the died individual 
to possibly colonize the empty site. In the networked landscape the neigh- 
borhood of a site is defined by the closest upstream and downstream sites. 
Limited dispersal promotes the clumping in space of species, which enhances 
their coexistence and survival probability [16, 33]. Indeed we find that in all 
the considered landscapes, persistence-time distributions still follow a power- 
law behavior characterized by smaller, non-trivial scaling exponents (namely 
a = 1.92 for the 3D, a = 1.82 for the 2D, a = 1.62 for the OCN, a = 1.50 
for the ID landscape, Fig. [3]) limited by an exponential cut-off. Remark- 
ably, the exponent obtained via simulation in a two-dimensional landscape 
(a = 1.82 ± 0.01) is close to those found in both breeding birds and herba- 
ceous plants datasets (a = 1.83 ± 0.02 and a = 1.78 ± 0.08, respectively). 

We also study how persistence-time distributions deducted from the the- 
oretical model change with dispersal broader than nearest neighbors (see 
Fig. S3). As expected, as long as the mean dispersal distance remains small 
with respect to the system size, the distribution eventually ends up scaling as 
the one predicted by the nearest-neighbors dispersal. We also relax the neu- 
tral assumption by implementing an individual-based competition/survival 
tradeoff model [16J. Specifically, species with higher mortality rates are as- 
sumed to hold less competitive ability to colonize empty sites [21 135] . It is im- 
portant to note that the trade-off model also exhibits power-law persistence- 
time distribution with exponents indeed close to those shown by the neutral 
model (see Fig. S4). Our theoretical results are thus robust with respect 
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both to changes in the dispersal range and to relaxations of the neutrality 
assumptions. This confirms our expectation that a power-law distribution for 
species persistence times is the result of emergent behaviors independent of 
fine ecological details, thus supporting the neutral assumption that effective 
interaction strength among species is weak [SB] and does not significantly 
constrain the dynamics of ecosystems. We also note that our results are 
not seen as a test for the neutrality hypothesis for breeding birds or herba- 
ceous plants dynamics, but rather as tools to reveal emerging universal and 
macroscopic patterns [37J EH] • 

4 Discussion and Conclusions 

In the previous section we have established a hierarchy of scaling exponents 
ranging from the smallest, proper to one-dimensional (1-D) interactions, to 
larger values namely for directional (network-like), 2-D, and 3-D dispersals. 
We thus suggest that the coherence of the empirical scalings would stem from 
the two-dimensional isotropic nature of the environmental matrix available 
to the ecological processes relevant to both breeding birds and herbaceous 
plants. 

We also suggest that species persistence-time distribution, owing to its 
robustness and scale-invariant character, is a synthetic descriptor of ecosys- 
tem dynamics and of biodiversity. In fact, other key macroecological patterns 
are intimately related to the persistence-time distribution. A first clear ex- 
ample is the direct link with ecosystem diversity, as explained below. In our 
framework species emerge as a point Poisson process with rate A = uN and 
last for a persistence time r. The mean number of species S in the system 
at a given time is therefore S = A(r) [39J where (r) is the mean persistence 
time. Therefore, the smaller exponents found, say, for networked environ- 
ments with respect to two-dimensional ones, imply longer mean persistence 
and, in turn, higher diversity. This echoes recent results suggesting a higher 
diversity of freshwater versus marine ray-finned fishes [1H1 SI] ■ 

Another evidence of the effective way in which species persistence-time 
distribution can characterize ecosystem diversity is the link with the species- 
area relationship, which characterizes the increase in the observed number 
of species with increasing sample area. The spatial extent of the breeding 
bird dataset and the tools developed for the data analysis allow us to study 
how the persistence-time distribution depends on the spatial scale of analysis 
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(Fig. |4]A) . As expected, while the scaling exponent remains the same, the 
diversification rate v decreases with the geographic area A and is found to 
closely follow a scaling relation of the type v oc , with (5 = 0.84 ± 
0.01 (Fig. [4j3) , for a wide range of areas. This scaling form of the cut-off 
timescale 1/u can be related to the species-area relationship. Assuming that 
the number of individuals scales isometrically with the sampled geographic 
area [8j [2], i.e. N oc A, and given that (r) = J tp T (t)dt oc z/*~ 2 (see SI) one 
gets: 

S = A (r) oc A 1 -^^ = A z . (2) 

The observational values /3 = 0.84±0.01 and a = 1.83±0.02 give an exponent 
z = 0.30 ± 0.02 which is close to the species-area relation measured directly 
on the data for the same range of areas (z = 0.31±0.02, Fig.|4p). Conversely, 
one could have used the observed species-area exponent to infer the scaling 
properties of the diversification rate. 

Finally, from a conservation perspective, a meaningful assessment of species' 
local extinction rates is deemed valuable. We propose the distribution of the 
times to local extinction r e (Fig. [2]A.) as a tool to quantify the dynami- 
cal evolution of the species assembly currently observed within a given geo- 
graphic area. Mathematically, r e is defined as the time to local extinction of 
a species randomly sampled from the system, regardless of its current abun- 
dance. When Eq. [T] holds for persistence times, the distribution of the times 
to local extinction p Te (t) is shown to scale as p Te (t) oc t x ~ a e~ vt (see Materials 
and Methods). Therefore, not only do the developed theoretical and opera- 
tional tools allow to infer the scaling behavior of persistence times, but also of 
the time to local extinction even from relatively short observational windows. 
Although these patterns cannot provide information about the behavior of a 
specific species or of a particular patch inside the ecosystem considered (e.g. 
a biodiversity hot-spot) they can effectively describe the overall dynamical 
evolution of the ecosystem diversity. In particular the scaling behavior al- 
lows to extrapolate species persistence-time distributions for wide geographic 
areas, which are hard to estimate, from measures of persistence on smaller 
areas, which are, on the contrary, more practical and feasible. We thus con- 
clude that the biogeographical characters of species persistence, stemming 
from the structure of the spatial interaction networks and from local con- 
straints to species emergence rates, add a new ingredient to a rich literature 
bearing major implications for the inventory of life on Earth. 
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5 Materials and Methods 



5.1 Inference of the Persistence-Times Distribution from 
a Finite Observational Period 

The exact derivation of the probability distribution of the variables r' and 
t" (Fig. follows. In this theoretical framework, The probability vdt 

of observing a diversification event in a time step dt is assumed to be a 
constant, thus species emergence in the system due to migration or speciation 
is seen as a uniform point Poisson process with rate A = vN (where N is 
total number of individuals in the system and A has the dimension of the 
inverse of a generation time). We term to the emergence time of a species 
in the system, and T and Tf = T + AT W the beginning and the end of the 
observational time window, respectively. A species emerged at time t will 
be continuously present in a geographic region for its persistence time r until 
its local extinction at time to + r. 

We first analyze the distribution of t", the most complex case. The 
variable t" can be expressed as function of the random variables r and to, 
which are probabilistically characterized. We can distinguish four different 
cases (Fig. [2}V): 

1. the species emerges and goes locally extinct within the time window; 

2. the species emerges during the observations and it is still present at the 
end of the time window; 

3. the species emerges before the beginning of the observations and goes 
locally extinct within the time window; 

4. the species is always present for all the duration of the observations, 
or, mathematically: 

{t, if T < to < Tf and t + t < Tf 

Tf - t , if T < t < Tf and t + r > T f 

to + r-To, if < t < T and T <to + T<T f 

T f - T , if < t < T and t + r > T 

We express the probability of observing r" conditional on a persistence 
time of duration r as: 
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p T "(t\r) = 

± ((6(t - t)Q(t - r o )0(T / - (t + r))0(T / - T - r)) + 
+ (S(Tf -t - t)e(t - T o )0(T / - t )e(t - (T f - r))) + 
+ (6(t + r - T - t)Q(t )e(T f -to- r)6(T - t )e(t - T + r)) + 

+ (S(T f -To- t)e(t )e(T - to)e(t - + r)e(r -T f + T ))) , (3) 



where the operator (•) is the ensemble average with respect to the random 
variable to, S(x) and Q(x) are the Dirac delta distribution and the Heaviside 
function, respectively. is the normalization constant. 

When comparing analytical and observational distributions, we assume 
that the system is at stationarity and unaffected by initial conditions, i.e. 
T Q is far from the beginning of the process. Mathematically this is obtained 
taking the limit To, Tf — > +00 with Tf — To = AT W . By solving the ensemble 
averages and by marginalizing with respect to r, Eq. [3] finally takes the form 
(see SI for a step by step derivation): 



Pr»(t) = ^ ( (AT W - t)p T {t)Q(AT w - t) + 



+ e{AT w -t) / p T (r)dr + 
Jt>o 

POO 

+ e(AT w -t) / p T (T)dr + 



t>0 



+ 5(t-AT w )[ (r - AT w )p T (r)dr) , (4) 
Jat w / 

where M simplifies to: 

N = AT W + (r) - 2 AT W P T (AT W ) + 2 (J (Pr(t) - tp T (t))dt) , 

with P T (t) = f^ 00 p T (r)dr being the exceedance cumulative distribution of 
the persistence-time probability density function. 
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The variable r' comprises only the first of the four cases listed in Eq. |3j 
Thus the probability distribution p T '(t) follows directly from the first term 
of Eq. H 

p T ,(t) = -^{AT w -t)p T (t)e(AT w -t), 
where the normalization constant M' is equal to 

W = / (AT w -r)p T (r)dr, 
Jo 

which completes the derivation. 



5.2 Distribution of Times to Local Extinction 

We term r e the time to local extinction of a species randomly sampled among 
the observed assembly at a certain time T (Fig. |2]A.). Analogously to the 
derivation described above, we can express r e as: 

Te = t + T-T if < t < T and t + r > T. 

We then express the probability distribution of the times to local extinc- 
tion conditioned to a persistence time r as: 



PtMt) = \(8{t - (f„ + r - r))e(f„ + t - T)Q(T - f )e(* )>, 

where the constant C ensures proper normalization. Solving the ensemble 
average operators yields 

p Te (t\r) = ^(r-T)Q(t-T + T). 

Marginalizing over r and considering the system at stationarity (T — > 
+00), we finally obtain 

1 [°° 

where C is simply (r). Eq. |5] allows to derive the distribution of the times 
to local extinction given the persistence-time distribution. Particularizing 
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now to the case of persistence-time distributions of the shape p T (t) oct a e ut , 
Eq. [5] translates into p Te (t) oc t l ~ a e~ vt . 
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Figure 1: Species persistence times. Persistence time r within a geo- 
graphic region is defined as the time incurred between a species' emergence 
and its local extinction. Recurrent colonizations of a species define different 
persistence times. The number of species in the ecosystem as a function of 
time (gray shaded area) crucially depends on species emergences and per- 
sistence times. We analyze two long-term datasets about North America 
breeding birds |22j and herbaceous plants from Kansas prairies [23]. The 
inset shows the observational routes of the Breeding Bird Survey. Aggregat- 
ing local information comprised in a given geographic area, we reconstruct 
species presence-absence time-series that allow the estimation of persistence- 
time distributions. 

Figure 2: Empirical persistence-time distributions, (a) A schematic rep- 
resentation of the variables that can be measured from empirical data over 
a time window AT W : r', persistence times that start and end inside the 
observational window, and r" , which comprises r' and all the portions of 
persistence times seen inside the time window that start or/and end outside. 
Times to local extinction r e are also presented, (b) Breeding birds and (c), 
herbaceous plants probability density function p(t) of r' (green), r" (blue) 
and persistence time r (red). Filled circles and solid lines show observational 
distributions and fits, respectively. The best fit is achieved with p T (t) oc t~ a 
with a = 1.83 ± 0.02 and a = 1.78 ± 0.08 for breeding birds and herbaceous 
plants, respectively. Note that previous estimates [6J for (b) are revisited 
here in the light of the new tools and of a longer dataset. The spatial scale 
of analysis is A = 10, 000 km 2 and AT W = 41 years for (b) and A=l m 2 and 
AT W = 38 years for (c). The finiteness of the time window imposes a cut-off 
to p T '(t) and an atom of probability in t = AT W to p T »(t), which corresponds 
to the fraction of species that are always present during the observational 
time. p T (t) and p T '{t) have been shifted in the log-log plot for clarity 

Figure 3: Persistence-time distributions are dependent on the structure 
of the spatial interaction networks, (a) Persistence-time exceedance proba- 
bilities P T (t) (probability that species' persistence times r be > t) for the 
neutral individual-based model [151 US] with nearest-neighbor dispersal im- 
plemented on the different topologies shown in the inset. Note that in the 
power-law regime if p T (t) scales as t~ a , P T (t) oc t~ a+1 . The scaling exponent 
a is equal to 1.50±0.01 for the one- dimensional lattice (red), a = 1.62±0.01 
for the networked landscape (yellow), 1.82±0.01 and 1.92±0.01 respectively 
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for the 2D (green) and 3D (blue) lattices. Errors are estimated through 
the standard bootstrap method. The persistence-time distribution for the 
mean field model (global dispersal) reproduces the exact value a = 2 (black 
curve). For all simulations v = 10 -5 and time is expressed in generation 
time units [TT]. Bottom panels sketch the color-coded spatial arrangements 
of species in a networked landscape (b), in a two-dimensional lattice with 
nearest neighbor dispersal (c), and with global dispersal (d). 

Figure 4: Biogeography of species persistence time, (a) Observational 
distributions p T >(t) and p T "(t) (interpolated solid circles) for the breeding 
bird dataset and corresponding fitted persistence-time distributions p T {t) oc 
t~ a e~ ut (solid lines) for different scales of analysis: Area A = 8.5 • 10 4 km 2 
(green), A = 3.4 • 10 5 km 2 (blue), A = 9.5 ■ 10 5 km 2 (red). v(A) provides 
the cut-off for the distribution, whose scaling exponent is unaffected by ge- 
ographic area. Note that the position of the cut-off of p T (t) is inferred from 
the estimate of the atom of probability of p T "{t) which is more sensitive 
to the scale of analysis; (b), Scaling of the diversification rate v with the 
geographic area v oc A -13 , (3 = 0.84 ± 0.01; (c), Empirical species-area re- 
lationship (SAR). The plot shows the mean number of species S found in 
moving squares of size A. We find S oc A z , z = 0.31 ± 0.02. Slope and 
confidence interval have been obtained averaging 41 SARs, one per year of 
observation. 
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